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Abstract 

We study the emission of scalars into the bulk from a higher- dimensional rotating 
black hole. We obtain an analytic solution to the field equation by employing matching 
techniques on expressions valid in the near-horizon and far-field regimes. Both analytic 
and numerical results for the absorption probability, in the low-energy and low-angular 
momentum limit, are derived and found to be in excellent agreement. We also compute 
the energy emission rate, and show that the brane-to-bulk ratio of the energy emission 
rates for scalar fields remains always larger than unity in the aforementioned regime. 



1 Introduction 



Theories of extra spatial dimensions [U [2] offer hope for the resolution of the hierarchy 
problem arising between the scale of gravity and that of other fundamental interactions. 
Within this framework gravity, and possibly scalar fields, propagate in the full (4 + Tri- 
dimensional continuum (bulk), while ordinary matter is trapped on a four-dimensional 
hypersurface (brane). The resolution of the hierarchy problem arises because the fun- 
damental scale of higher- dimensional gravity can be close to the other particle physics 
scales, while the traditional scale of gravity (Planck mass) is related to it in terms of the 
volume and number of extra dimensions. In this context, trans-planckian collisions can 
give rise to potentially observable black hole creation [3]. These higher dimensional black 
holes [4J could be produced either in ground-based colliders or in cosmic ray interactions 
and detected through the emitted Hawking radiation [5] . There is a considerable amount 
of recent literature on the subject (see [6]) including both numerical and analytic studies 
of the Hawking radiation in the Schwarzschild phase [H El El HOl [HI [121 [131 HH [151 US] 
as well as the spin-down phase [I71[l8lQIl[IB[IDl[22l^ 

In previous articles we provided analytic results for the evaporation of higher- 
dimensional black holes through the emission of scalars [27] , fermions and gauge bosons 
[28] on the brane. Here, we continue our program of analytic study by focusing on the 
emission of scalar degrees of freedom in the bulk from rotating higher- dimensional black 
holes. In section 2 we set up the general theoretical framework, considering the gravita- 
tional background corresponding to a rotating higher- dimensional black hole and writing 
down all relevant equations for a scalar field in such a spacetime. In section 3 we pro- 
vide an analytic solution of the radial equation using a well-known matching technique. 
This technique consists of, first, deriving the solution in the near horizon regime, then 
deriving the equivalent far field limit, before finally stretching and matching the two 
forms in an intermediate zone. In this way an analytic expression for the radial part 
of the field valid throughout the entire spacetime is constructed. The solution obtained 
in this way is valid in the low-energy and low-angular-momentum approximation. Our 
solution is used in section 4 to calculate the absorption probability which characterises 
the emitted Hawking radiation. We also perform a complimentary numerical analysis, 
valid in the same regime, and compare the two sets of results. Plots are depicted for 
various dimensional spacetimes and angular momentum modes. In section 5 we compute 
the energy emission rate and the brane-to-bulk ratio of the energy emission rates, and 
provide corresponding plots. Finally, in section 6 we present our conclusions. 



2 General Framework 

The gravitational background around a (4+n)-dimensional, neutral, rotating black hole is 
described by the well-known Myers- Perry solution [4J. Here we will be interested in black 
holes created by an on-brane collision of particles, and may therefore assume that the 
corresponding metric will depend only on one non-zero angular momentum component, 
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and dVt\ (9i, 92, ■ ■ ■ , # n _i, 0) is the line-element on a unit n-sphere. The black hole's mass 
Mbh and angular momentum J are related to the parameters a and fi as follows 
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with G being the (4 + n) -dimensional Newton's constant, and A n+2 the area of a (n + 2)- 
dimensional unit sphere given by A n+2 = 27r(" +3 )/ 2 /r[(n + 3)/2]. Finally, the black hole's 
horizon radius follows from the equation A(r^) = 0, and is found to satisfy the 
equation: r^ +1 = /i/(l + a 2 ), where a* = a/r^. 

In this work we will study the emission of Hawking radiation, in the form of scalar 
fields, in the higher-dimensional space. We therefore need to consider the equation of 
motion of a massless scalar field propagating in the gravitational background ([T]). This 
is given by 
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The above equation can be separated by assuming the factorised ansatz 

$ = V m ^(r) 5(0) F jn (^, . . . , n _ 1; 0) , 
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where Yj n {9\, . . . , n _i, 0) are the hyper spherical harmonics on the n-sphere that satisfy 
the equation [29] 
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The functions i?(r) and S(9) are then found as the solutions to the decoupled equations 
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respectively, that first appeared in the literature in Ref. [3U]. In the above, 

K = (r + a ) oj — am , Ajg m = Ej£ m + a uo — 2amu . (10) 

The angular eigenvalue E^ m {aio) provides a link between the angular and radial 
equation and, as in the case of on-brane emission [27J OS EH EH El] , there is no closed 
analytic form for its value. It may however be expressed as a power series in au [3T] 132] . 
For the purpose of our analysis, valid only in the low-u and low-a limit, it suffices to keep 
a finite number of terms and we will therefore truncate the series at the 5th ordeJE In 
order for the power series to converge in the limit au —>■ by terminating at finite order, 
a number of restrictions are imposed on the allowed values of the integer parameters 
(j, £, m) specifying the emission mode : in general, m may take any integer value and j 
and I any positive or zero integer value providing [31] 

£>j + \m\ and - - (j + |m|) £ {0,Z+}. (11) 

By using the power series form of the angular eigenvalues, we may now proceed to 
solve Eq. ([8]) analytically. The solution for the radial function R(r) will help us determine 
the absorption probability |v4^ m | 2 for the propagation of a massless scalar field in the 
bulk, a quantity that characterises the Hawking radiation emission rates of the black 
hole. 



3 Analytic Solution 



In order to obtain a solution to the radial equation ([8]), we will use a well-known ap- 
proximate method: we will first solve the equation close to the horizon of the black hole 
(r ~ r/i), and then far away from it (r ^> r^). Finally, we will smoothly match the two 
solutions in the intermediate zone, thus creating an analytical solution for the whole 
radial regime. 

Starting from the near-horizon regime, we perform the following transformation of 
the radial variable 127, 281 
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where A(r) = {n + 1) + {n — 1) a 2 /r 2 . Then, near the horizon (r ~ r^), Eq. ([8]) takes the 
form 
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1 The expression we use for Ej£ m is based on the analysis of [31] but disagrees slightly with the 
version given there as the sign of the second order term is reversed. This is necessary so that, in the 
limit j,n — > 0, the expression for Ej£ m correctly reproduces the on-brane eigenvalues that have appeared 
in the literature previously. 
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where A* and K* are given by 

A* = (n + 1) + (n- l)a 2 , if* = (1 + a*)a;* - a*m , (14) 

and 

4a 2 

A^l-^. (15) 

By employing the transformation 

JW/) = A-r (1 - /)' ^(o, b, c; f) , (16) 

equation ( fl3l) can be brought to the form of a hypergeometric differential equation [33J, 
with a = a + (3 + D* — 1, 6 = a + /5, and c = 1 + 2a. In addition, the parameters a and 
/3 are given by 

« = ±— (17) 
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The general solution of the radial equation ffl3|) may be written in terms of the hyper- 
geometric function F as follows 

R NH U) = A„ r(l-ffF(a,b,c; f) 

+ A + f- a (l-f)PF(a-c + l,b-c + l,2-c;f). (19) 

In this general near-horizon solution, we must now impose the boundary condition that 
no outgoing modes exist near the black hole's horizon. In the limit r — > r% we get 
/(r) — > 0, and the near-horizon solution (fl9|) becomes 

ifon (/) - A_ / ±iX */ A * + /=f*^M« = A_ e ±ifcy + e Ti ^ , (20) 

where 



k = u) — mVt — u 5 r , (21) 

+ a z 



and |/ is a tortoise-like coordinate defined by y — r h {l + a 2 ) ln(/)/A*. By imposing this 
boundary condition, we can set either A_ = or A + = 0, depending on the choice for 
a. The two choices are clearly equivalent, so we choose a = a_, which imposes A + = 0. 
This brings the near-horizon solution to the final form 

R NH (f) = A. f a (1 - ff F(a, b, c; f) . (22) 

Finally, the convergence criterion for the hypergeometric function, i.e. Re(c — a — b) > 0, 
should also be applied, leading us to choose /? = /?_. 

We now focus our attention on the far-field regime, r > r/,. In this limit, the sub- 
stitution R(r) = r~(^ _ )i?(r) brings Eq. (jSJ) into the form of a Bessel equation [33], in 
terms of z = our, 
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Thus, the solution in the far-field regime may be written as 

R FF (r) = 4k J » M + -Sr Y » M > ( 24 ) 

r 2 t 2 

with J v and Bessel functions of the first and second kind, respectively, and v = 

In order to construct a full analytic solution valid for all r, we must smoothly match 
the two asymptotic solutions (1221) and (1211) in an intermediate regime. To this end, we first 
focus on the near- horizon solution (122)) . and shift the argument of the hypergeometric 
function from / to 1 — / by using the following relation [33J 
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The function f(r) may be alternatively written as 
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In the limit r 3> r^, the (r/r^) 2 component in the denominator of the second term 
dominates, and fir) goes to unity for n > 0. Then, the near- horizon solution (125)) can 
be written as 

i^(r) ~ A 1 r- ( - n+ W + A 2 r( n+ W +D *-y , (27) 
with Ai and A 2 defined as 

Al " A -K 1 + «J^ J r( c -a)r( c -6)' 
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Next, we expand the far- field solution to small values of r, by taking the limit r — ► 
in Eq. flM} : 

J2 M~^iM- * 2 rM (29) 
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Then, taking the small a* and cu* limit in the power coefficients of r - so that we can 
ignore terms of order a 2 , a*a;*) or higher - we can achieve exact matching since 
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We would like to stress here that, in order to achieve a higher level of accuracy in our 
analysis, no expansion is performed in the arguments of the gamma functions, and terms 
to order (acu) 4 are retained in the expansion of the eigenvalues. Then, the matching of 
the two asymptotic solutions leads to the constraint 
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The above guarantees the existence of a smooth, analytic solution for the radial part of 
the wavefunction for all r, valid in the low-energy and low-rotation limit. 



4 The Absorption Probability 



With the solution to the radial equation (jSJ) at our disposal, we may now compute the 
absorption probability. To this end, we expand the far- field solution fl23|) for r — > oo. By 
using standard formulae for the Bessel functions [33J, we find 
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We see that for large distances from the black hole, the solution reduces to an incoming 
and an outgoing spherical wave, allowing us to compute the absorption probability from 
the ratio of their amplitudes: 
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The above equation, in conjunction with Eq. fl3Tl) . is our main analytic result for the 
absorption probability characterising the emission of massless scalar fields in the bulk, 
from a rotating, uncharged black hole, in the low-energy and low-angular momentum 
limit. 

In Fig. ([T]) we plot the absorption probability for the first partial waves for n = 2, 
a* = 0.4, with the values of j, £, m obeying the restrictions fill I) . One may easily observe 
the dominance of the first partial wave j = i = m = over all others, and the suppression 
of \Aje m \ 2 as the values of the angular momentum numbers increase. In the plot, we also 
see the appearance of superradiance [M] for modes with positive m, where the absorption 
probability takes negative values. Figure (jTJ) actually depicts two sets of curves: the first, 
denoted by solid lines, gives the value of \Aj£ m \ 2 that follows from our analytic result (J3"3"j) . 
while the second set, denoted by dashed lines, gives the value following from integrating 
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Figure 1: Absorption probabilities |«4j\e m | f° r a bulk scalar field, for n = 2 , a* = 0.4, 
and various combinations of (j,£,m). 



the radial equation ([8]) numericalljo As in the case of scalar [27] and higher-spin fields 
[2"g] propagating on the brane, our approximate analytic method leads to results that are 
in excellent agreement with the exact numerical ones, not only in the low-energy regime 
but beyond this also. 

Focusing on the dominant first partial wave, in Fig. [2](a) we demonstrate the depen- 
dence of the absorption probability on the rotation parameter a*, for fixed n — 5. One 
may clearly see that an increase in the rotation of the black hole causes an enhancement 
in the value of | A,-£ m | 2 in the low- and intermediate-energy regimes. In Fig.[2](b), we inves- 
tigate instead the effect of the number of extra dimensions on the value of the absorption 
probability, while we fix a* = 0.5. As is clear from the plot, the value of \Aje m \ 2 in the 
low-energy regime is strongly suppressed as n increases. The same behaviour in terms of 
n was found for bulk scalar fields propagating in the background of a higher-dimensional, 
spherically-symmetric black hole [8]. 

An interesting question is how the absorption probabilities for brane and bulk scalar 
fields in a rotating black-hole background compare. By examining the results presented 
in this section and in [27], one may easily conclude that the absorption probabilities for 
both species of scalar fields are enhanced as the black hole rotation parameter increases. 
In contrast to this, the value of the absorption probability increases with n for brane 
scalars, while it decreases for bulk scalar fields. Important conclusions can also be drawn 
by directly comparing Fig. [T] with the corresponding figure in [27J: for the same values 
of n and a*, the absorption probability for brane scalar fields is consistently larger than 
that for bulk scalars by almost 3 orders of magnitude, both for superradiant and non- 
superradiant modes. The same observation was made in (23] in the 5-dimensional case - 
here, we have shown that this behaviour persists for higher values of n also. 

2 In the numerical integration of Eq. ([5]) the power series expansion of the eigenvalue Ejg m |31j up to 
5th order, valid in the low-aw limit, was again used. 
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Figure 2: Absorption probabilities for the bulk scalar mode j 
n = 5 and various a*, and (b) a* = 0.5 and various n. 
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A compact analytic expression for the absorption probability may be derived in the 
very low-energy limit. For uj — > 0, we obtain B oc l/u +n+1 , therefore 
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Substituting for B from Eq. (1311 . we take 
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Focusing our attention on the dominant mode j = i = m = 0, and performing an 
analysis similar to that in [27], where the absorption probability of the dominant scalar 
mode on the brane was found in the same limit, we obtain the result 
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The above result allows us also to compute the absorption cross-section o"o for the 
dominant scalar bulk mode in the asymptotic low-energy regime. By adapting the formula 
of Ref. [35] to the geometrical set-up of our analysis, we may write 
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Figure 3: Absorption cross-section for the bulk scalar mode j = I = m = 0, for n = 2 
and various values of a*. 

are the multiplicity of the jth partial wave in the expansion of the wave function over the 
hyperspherical harmonics on the n-sphere, and the horizon area of the (4+n)-dimensional 
rotating black hole, respectively Substituting for the absorption coefficient, we obtain 



For a* — ► 0, the above reduces to Ah, a behaviour that was found in [8]. For a* ^ 0, the 
numerical results (dashed lines) presented in Fig. [3] confirm that, also here, the low-energy 
limit of the cross-section tends to the area of the corresponding rotating black hole. The 
solid lines demonstrate the agreement of our analytic results with the numerical ones for 
small values of a*, and the expected deviation for large values of the rotation parameter. 

5 Energy Emission Rate 

Having found the absorption probability (Eq. I3"3"|) . we can now proceed to compute the 
rate of energy emission of massless scalar fields in the bulk. This is given by the expression 



The above differs from the 4-dimensional [36] expression in the presence of an additional 
sum over the new angular momentum number j, and from the 5-dimensional [T7] one in 
the introduction of the multiplicity of states Nj (|37|) following from the expansion of the 
wavefunction of the field in the n-dimensional sphere. The parameter k is defined in Eq. 
(1211) . while the temperature of the higher-dimensional, rotating black hole is 






(39) 




(40) 



9 



A useful check that may convince us of the validity of the above emission rate is 
to take the non- rotating limit. Then, Eq. (|3T?|) for the energy rate should reduce to the 
well-known result for bulk scalar emission from a (4+n)-dimensional Schwarzschild black 
hole [7J [8] 

d 2 E 1 v-^ u „ T . . l2 

vr = r E r it 1 i ^ ^ ' 41 

atrfu; 27T ' exp [w/T H J — 1 

where now Ni is the degeneracy of the £th mode of the harmonics on the (n + 2)-sphere 

(K + n + l)(/ + n)l 
£ £!(n + l)! 1 J 

In the limit a — > 0, we get k = u, and the absorption probability becomes independent 
of both m and j, retaining a dependence only on the principal quantum number I. Then, 
Eq. ( 139]) matches Eq. fj4Tl) providing the following relation holds 

In the second part of the above equation, we have used the fact that, according to the 
restrictions (fTTj) imposed on the quantum numbers, for each value of (j,£), m may take 
I — j + 1 values, and for each value of £, j may take the values < j < i. In order to 
prove Eq. (H3|) , we rewrite the factor 2j + n — 1 as (j + n — 1) + j , and split the sum in 
two parts. Then, if we further replace the index j by % — n + 1 in the first sum and by 
% — n + 2 in the second, the middle part of Eq. fj4"3l takes the form 

e+n-i e+n-2 

(t + n-l-i + l)C(i,n-l)+ ^ (£ + n - 2 - i + 1) C(i,n - 1) , (44) 

i=n— 1 i=n—l 

where C(s,r) is the combination function, C(s,r) = ^^yy. By using the identity 

s 

^2(s-i + l)C(i,r) = C(s + 2,r + 2), (45) 

the first sum in Eq. (j4"4"j) reduces to C(£ + n-|- 1, n+ 1), and the second to C(£ + n, n+ 1) 
- their sum can be easily shown to equal Ng. 

By using Eq. (|39l) . in Figs. H](a,b), we plot the energy emission rate for a higher- 
dimensional rotating black hole in the form of bulk scalar fields, as a function of the 
energy parameter wr/,, and in terms of the angular momentum parameter and number of 
extra dimensions, respectively. The profile exhibited by the absorption probability is also 
observed here: the emission rate is enhanced with a* in the low-energy regime, in analogy 
with results for brane scalars, fermions and gauge bosons [2QI (2U [2H [25], [26J [271 128] . 
while is suppressed in terms of n. Drawing experience from previous studies, we expect 
the enhancement in terms of a* to persist over the whole energy regime - on the other 
hand, it is quite likely, given the similarity of the results with those for bulk scalar fields 
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Figure 4: Energy emission rates for bulk scalar fields (a) for n = 5 and various a*, and 
(b) for a = 0.5 and various n. 



in a non-rotating background [8] , that the low-energy suppression with n will be replaced 
by a strong enhancement at the high-energy regime. 

We should finally address the question of the brane-to-bulk ratio of the rate of scalar 
field energy emission from a higher- dimensional rotating black hole. The answer to this 
question is important as it will define the amount of energy spent by the black hole in 
the observable brane channel. Whereas the energy emission rate for bulk scalars is given 
by Eq. (l39l) . the corresponding one for brane scalar degrees of freedom is [18j E3 EfE [21] 

d 2 E 1 UJ 2 

dtd^ = 2^f^ exp[Jfc/r H ]-l 1 eml ' ( } 

The comparison between the bulk and brane absorption probabilities discussed in the 
previous section has given a clear signal as to which emission is dominant, however the 
final comparison should involve the total emission rates where the different multiplicities 
of states have been taken into account. By using Eqs. fl39l) and (|4*6|) . the brane-to-bulk 
ratio for scalar fields from a rotating black hole is depicted in Figs. [5](a,b) in terms of the 
parameters a* and n. We may easily observe that although the exact value is both a* and 
n-dependent, the ratio of the brane to bulk emission rates always remains above unity, 
rendering the brane channel dominant. We remind the reader that a similar conclusion 
was also drawn in the case of a Schwarzschild-like higher-dimensional black hole [8]. 



6 Conclusions 

In this work we have performed an analytic study of the emission of scalar fields in the 
bulk from a higher- dimensional rotating black hole. By solving analytically the radial 
part of the equation of the scalar field, we have constructed a smooth solution, valid 
in the low-energy and low-angular-momentum regime. We were then able to compute 
the corresponding absorption probability, and to examine its behaviour in terms of the 
particle's angular momentum numbers and spacetime properties. The lowest scalar mode 
was found to be dominant, with its absorption cross-section equal to the horizon of the 
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Figure 5: Brane-to-bulk ratio of the energy emission rates for scalar fields (a) for n — 1 
and various a*, and (b) for a = 0.5 and various n 



higher-dimensional rotating black hole in the extreme low-energy limit. We demonstrated 
that the absorption probability for bulk scalar fields is enhanced with increasing angular 
momentum of the black hole and suppressed by the number of extra dimensions. We also 
performed a numerical analysis, valid again in the low-cu and low-a* regime, and showed 
that our analytic results are in excellent agreement with the exact numerical ones even 
up to intermediate-energy regimes. 

We then computed the energy emission rate in the bulk, and showed that it exhibits 
the same behaviour as the absorption probability in terms of a* and n. Finally, we 
calculated the brane-to-bulk ratio of energy emission rates for scalar fields in a rotating, 
higher- dimensional background and found that, in the low-u; and low-a* regime, it is 
always larger than unity. This is due to the fact that the larger multiplicity of scalar 
states in the bulk cannot counteract the fact that the bulk absorption probability is 
approximately three orders of magnitude smaller than the corresponding brane value in 
the low energy limit. It is this dominance of the brane absorption probability that largely 
determines the preference of the black hole to emit scalar field energy on the brane. The 
complete radiation spectrum - that may follow only from an exact, numerical analysis 
involving both the angular and radial parts of the field equation - is necessary to decide 
whether the brane-to-bulk ratio remains above unity over the entire frequency range. 
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